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Overview 


• MGA-1DSM 

• Developed Algorithms 

• Test Problems 

• Results 


MGA-1DSM 


• Multiple Gravity Assist with 1 Deep Space Maneuver 
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Phase-Wise Gene Recombination 


Decision Variables 
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Solution A 
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New Trial Solution 



Phase Genetic Solver 


• Single trial evaluations 


Trial creation 


• Trial creation by Phase-wise GA-style 
or DE-inspired recombination 

• "Bin" repository structure - requires 
an initialization period 
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Locally optimize 


• Non-exclusionary ''Kill Distance" 

• Population collapse mechanic 

• Main loop 

Creation - Probabilistic switch between 
GA and DE creation types 

Locally optimize 

Submit to repository 

Repeat 






Fitness 



Phase Genetic Solver - Population 
Collapse 


solutions 

• Fitness avalanche observed if necessary elements exist in population 

• If such material is missing, the bin convergences in fitness space 

• Proliferation of breakthrough genetic information moves evolution out of the 
local optimum 
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Test Problems 


• ESA GTOP Database inspired problems 

• Cassini 2 

• Messenger Reduced 

• Real-world problems 

OSIRIS-REx with 90° spacecraft target Sun angle 
OSIRIS-REx with 45° spacecraft target Sun angle 


Cassini 2 


• Academic model problem based on the Cassini mission 

• Bounded time of flight 

• Single journey from Earth to Saturn with a fly-by sequence of: 
Venus, Venus, Earth, Jupiter 

• Rendezvous condition at Saturn 

• Objective function: minimize Av 


Results - Cassini 2 


• Total Av : 8.4878 

• Flight time (y): 8.0000 
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Event # 6: 
unpowered flyby 
Earth 
8/17/1999 
r x = 15.387 km/s 
DEC = -15.0° 
altitude = 303 km 
m = 1990 kg 


Event # 1: 
launch 
Earth 
5/7/1997 
C3 = 0.053 km' /s 2 
DLA = -21.2 3 
m = 6100 kg 


Event # 2: 
chemical burn 
deep-space 
10/31/1997 
Ar = 3.069 km/s 
m = 2294 kg 
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Event # 5: 
unpowered flyby 
Venus 
6/18/1999 
r x = 9.364 km/s 
DEC = -1.8 
altitude = 3529 km 
m = 1990 kg 


Event # 7: 
unpowered flyby 
Jupiter 
1/30/2001 
r x = 9.724 km/s 
DEC = -3.7° 
altitude = 8048441 km 
m = 1990 kg 


0.0 


Event # 3: 
unpowered flyby 
Venus 
4/19/1998 

v x = 6.021 km/s 

DEC = 16.4 
altitude = 303 km 
m = 2294 kg 


Event # 8: 

rendezvous 

Saturn 

5/7/2005 

Ar = 4.742 km/s 

m = 439 kg 


x (km) 


Event # 4: 
chemical burn 
deep-space 
12/11/1998 
Ar = 0.447 km/s 
m = 1990 kg 
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Tested Tunings - Cassini 2 


Test Name 

Minimum Fitness 

Maximum Fitness 

Mean Fitness 

Standard Deviation of Fitness 

NM max iterations = 1000 
NM max iterations = 250 & 

8.507702203 

13.80813145 

10.19337892 

2.381004889 

NM stagnation limit = 30 

8.495142971 

13.8081 1605 

9.129719809 

1.583863494 

Disabled proximity kill 

8.487794692 

15.4574581 

9.517611548 

2.019381623 

NM initial step = 0.1 

8.519776442 

13.80824377 

12.40553514 

2.155671733 

NM repeat = 3 

8.491886681 

13.8081 1288 

9.534945467 

2.055517186 

NM initial step = 0.01 

8.529650524 

13.80819269 

9.742230666 

1.56202128 

Kill distance = 1 .25 

8.548351326 

13.80855295 

9.618825493 

2.094259801 

NM expansion parameter = 5 

8.49321 1225 

13.80810881 

9.213243008 

1.572387353 

Kill distance = 0.45 

8.489449107 

15.64071313 

10.84552859 

2.550076981 

Kill distance = 1.5 

8.548151998 

8.6373963 

8.588254452 

0.028344059 

DE creation percentage = 60 

8.510255988 

13.80812768 

10.37414168 

2.37686051 

NM max iterations = 250 

8.502856018 

13.80823687 

10.28526561 

2.330492865 

Fitness penalty weight = 50 

8.510825618 

9.445452757 

8.696852204 

0.354515305 

DE creation percentage = 30 

8.943567626 

16.04331072 

13.08227751 

2.115409215 

DE creation percentage = 80 

8.518332019 

8.578794109 

8.541124301 

0.015297718 

NM stagnation limit = 50 

8.500648814 

13.80810493 

9.037603714 

1.590176354 

Fitness penalty weight = 10 

8.50422969 

44.25559316 

17.5745441 1 

11.51069538 

Default 

8.501509161 

15.64214162 

1 1.20484013 

2.95955532 

DE creation percentage = 70 

8.517078919 

8.556319048 

8.534003909 

0.009798648 

Kill distance = 2.0 

8.638286021 

13.87964676 

10.2763394 

2.346441835 

Pareto mutation percentage = 100 

8.501850121 

11.16184864 

8.972233049 

0.816102361 


Messenger Reduced 


• Academic model problem based on the Messenger mission 

• Bounded time of flight 

• Single journey from Earth to Mercury with a fly-by sequence 
of: Earth, Venus, Venus 

• Rendezvous condition at Mercury 

• Objective function: minimize Av 



Results - Messenger Reduced 


• Total Av : 7.8483 


• Flight time (y): 4.1872 
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Results - Messenger Reduced 


• Total Av : 7.7005 

• Flight time (y): 
4.3597 


Event # 4: 
unpowered flyby 
Venus 
1/18/2000 
V x = 8.858 km/ a 
DEC = 12.3“ 
altitude = 2336 km 
m = 2159 kg 



Event # 5: 


unpowered flyby 
Venus 

' 

11/22/2001 

\ 

Hjo = 8.855 km/ s 

\ 

DEC = 2.5 = 

\ 

altitude = 7117 km 

\ 

m = 2159 kg 



Event # 1: 

launch 

Earth 

5/12/1998 

C3 = 1.180 km 2 /s 2 

DLA = 15.2“ 

m = 2959 kg 


Event # 3: 
unpowered flyby 
Earth 
8/3/1999 
r x = 4.108 km/ a 
DEC = 16.6“ 
altitude = 4900 km 
m = 2159 kg 


Tested Tunings - Messenger Reduced 


Test Name 

Minimum Fitness 

Maximum Fitness 

Mean Fitness 

Standard Deviation of Fitness 

DE creation percentage = 30 
NM max iterations = 1000 & 

7.700553833 

10.077949 

8.69836204 

0.957901076 

NM stagnation limit = 100 & 
Pareto mutation percentage = 100 
NM max iterations = 1000 & 

7.885601368 

9.838366048 

8. 1 3528245 1 

0.580954014 

NM stagnation limit = 100 

7.848294059 

9.196548589 

8.092378624 

0.424963406 

Disabled proximity kill 

7.883441903 

1 1.46388641 

9.370414423 

1.286407838 

NM max iterations = 1000 

7.885493493 

11.46414506 

8.830137907 

1.18791056 

NM repeat = 3 

7.884693973 

10.5763399 

8.818799601 

0.8565251 1 

DE creation percentage = 70 

7.891120342 

9.968710382 

8.988314296 

0.931023514 

NM stagnation limit = 50 

7.884667061 

9.950609605 

8.665714852 

0.95438583 

NM stagnation limit = 100 

7.884242668 

9.951 1 19188 

8.163303062 

0.632826718 

NM initial step = 0.01 

7.883914168 

1 1.46634384 

9.96979681 

1.206551117 


OSIRIS-REx 


• 2 Journey Mission 

Journey 1: Earth to Bennu with 1 Earth fly-by 
Journey 2: Bennu to Earth 

• Bounded allowable spacecraft target Sun angle at Bennu 

• Trial 1: 90° 

• Trial 2: 45° 

• Objective function: minimize Av 


Results - OSIRIS-REx 90° 


• Total Av : 0.7238 

Outbound Av: 0.4624 

• Flight time (y): 7.0041 Launch day: 9-23-2016 



1.0 - 


0.5 - 


5000 10000 15000 20000 25000 


[Tfest2 




Tested Tunings - OSIRIS-REx 90° 


Test Name 

Minimum Fitness 

Maximum Fitness 

Mean Fitness 

Standard Deviation of Fitness 

Disabled proximity kill 

0.723849797 

2.932235497 

1.173036182 

0.68831895 

NM initial step = 0. 1 

0.725279214 

1.796980154 

0.835174524 

0.320672698 

NM stagnation limit = 50 

0.724487126 

1 .233349354 

0.810434537 

0. 154609675 

NM max iterations = 1 000 

0.724595703 

1.608428143 

0.954556275 

0.291337706 

NM stagnation limit = 100 

0.724276918 

1.607709109 

0.852585628 

0.260491888 

Kill distance = 1 .2 

0.725769679 

1.288751624 

0.802798213 

0.16791975 

Kill distance = 0.45 

0.724017532 

1.607853071 

0.898080663 

0.258426999 

DE creation percentage = 30 
NM max iterations = 1 000 & 

0.72518409 

17.29093509 

3.433378663 

4.746951776 

NM stagnation limit = 100 

0.724384127 

873.0522079 

122.7224839 

257.1555035 

DE creation percentage = 70 

0.725050606 

1 .609679952 

0.831758675 

0.26446865 

Pareto mutation percentage = 100 

0.724466525 

1 .609754007 

0.869859487 

0.256750254 

Default 

0.725657145 

3.738515419 

1.107804668 

0.881684002 


Results - OSIRIS-REx 45° 


• Total Av : 0.7241 

Outbound Av: 0.4625 

• Flight time (y): 7.0040 Launch day: 9-23-2016 
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Event # 3: 

chemical burn 

deep-space 

4/26/2018 

Ar — 0.055 km/s 

m = 1908 kg 



Event # 6: 

chemical burn 

deep-space 

3/6/2023 

Ar — 0.012 km/s 

m = 1420 kg 



Tested Tunings - OSIRIS-REx 45° 


Test Name 

Minimum Fitness 

Maximum Fitness 

Mean Fitness 

Standard Deviation of Fitness 

Disabled proximity kill 

0.794681232 

3.797895363 

1 .64365942 

1.101750074 

Pareto mutation percentage = 100 

0.799880971 

6.160661388 

1 .468606482 

1.583540108 

NM max iterations = 1 000 

0.724436189 

4.045370978 

1.159974362 

0.987205563 

NM stagnation limit = 100 

0.724134498 

1 .06855567 1 

0.801919046 

0.124221252 

Default 

0.762818765 

1.64313778 

0.957427469 

0.248590476 

Kill distance = 1.2 

0.800867037 

1.834345561 

0.943784705 

0.309631387 

Kill distance = 0.45 

0.754246276 

5.67378172 

1.451568898 

1.443897955 

NM initial step = 0. 1 

0.79777825 1 

2.361 178032 

0.993787507 

0.461541481 

DE creation percentage = 70 
NM max iterations = 1 000 & 

0.798337374 

1.084989391 

0.879957468 

0.105577135 

NM stagnation limit = 100 

0.795939001 

0.96675381 

0.836301535 

0.065568903 

DE creation percentage = 30 

0.810500627 

13.98726773 

2.472446933 

3.860479861 

NM stagnation limit = 50 

0.80221576 

2.052751678 

0.966982707 

0.368824685 


Conclusions - FIX 


• The Phase Genetic Solver works very well on the chemical 
mission problems tested 

• Tuned problems are solved reliably 

• A tuning law has yet to be developed 

• This method can be extended to any problem where grouped 
decision variables have additional meaning 


Questions? 


